Bose-Einstein Condensation of Magnons in Cs 2 CuCl4 
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We report on results of specific heat measurements on single crystals of the frustrated quasi-2D 
spin-1/2 antiferromagnet CS2CUCI4 (Tn = 0.595 K) in external magnetic fields B < 12 T and for 
temperatures T > 30 mK. Decreasing B from high fields leads to the closure of the field-induced 
gap in the magnon spectrum at a critical field B c ~ 8.51 T and a magnetic phase transition is 
clearly seen below B c . In the vicinity to B c , the phase transition boundary is well described by the 
power-law T C {B) oc (B c — B) 1 ^ with the measured critical exponent <f> ~ 1.5. These findings are 
interpreted as a Bose-Einstein condensation of magnons. 
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In a quantum antiferromagnet (AFM) a fully spin- 
polarized state can be reached at high magnetic field B 
exceeding a saturation field B c . In this state, spin ex- 
citations are gapped ferromagnetic magnons. With de- 
creasing B and passing through B Cl an antiferromagnetic 
long-range order of the transverse spin component devel- 
ops. Provided the symmetry of the spin Hamiltonian 
is such that the rotational invariance around the applied 
field is preserved, the transverse spin component ordering 
can be regarded as a Bose-Einstein condensation (BEC) 
in a dilute gas of magnons. This concept was formu- 
lated theoretically many years ago Q, @ . For most of 
the known AFMs, B c can be well above 100 T. An ex- 
ceptionally low and easily accessible saturation field of 
B c w 8.5 T, however, is needed in the quantum spin- 
1/2 AFM CS2CUCI4. In this system the dominant ex- 
change spin coupling J is rather weak, J = 4.34(6) K 
The other isotropic spin coupling constants and the 
anisotropic Dzyaloshinsky-Moriya (DM) interaction are 
smaller and were determined with high accuracy by neu- 
tron experiments Q. Thus, the spin Hamiltonian in- 
volves the isotropic exchange Hq, the DM anisotropic 
term Hjjm and the Zeeman energy Hb and is given by 
H = Hq + Hdm + Hb- 

CS2CUCI4 falls into the class of easy-plane AFMs with 
?7(l)-rotational invariance around the crystallographic a- 
axis. Thus, for B applied along the a-axis, the U(l) sym- 
metry can be broken spontaneously due to the transverse 
spin component ordering at T c . This is accompanied by 
the appearance of a Goldstone mode with linear disper- 
sion, which is interpreted as signature of a magnon BEC 
P|. However, an unambiguous evidence for a BEC de- 
scription of the field-induced phase transition would be 
the determination of the critical exponent <j) in the field 
dependence of the critical temperature 



T C {B) cx (B c - Bf/ 



(1) 



Theory for a 3D Bose gas predicts a universal value 
0bec = 3/2 0, which coincides with the result of a 
mean- field treatment . 



A magnon BEC in TICUCI3 was recently reported 
0, Q. In this quantum AFM with a dimerized spin- 
liquid ground state, the saturation field is rather high, 
B C 2 ~ 60 T, and the BEC transition was studied near 
the first critical field, B cA ~ 5.6 T. At B = B cA , the 
singlet-triplet excitation gap is expected to close and a 
BEC occurs for B > B c \ [||. However, a few experimen- 
tal findings show deviations from a pure magnon BEC: 
An anisotropic spin coupling (of unknown nature) might 
produce a small but finite spin gap in the ordered state 
for B > B Ct \ [8j and the reported critical exponent <f> is 
somewhat larger than predicted by theory [6j. 

In this Letter we report on specific heat measurements 
9] on single crystals of CS2CUCI4 at low temperatures 
(30 mK < T < 6 K) and high magnetic fields (B < 
12 T) applied along the crystallographic a-axis. The aim 
of this thermodynamic study was (i) to trace the field 
dependence of T C {B) near B c , i.e., to extract the power 
law according to eq. ^ and (ii) to determine the closure 
of the spin gap. The access to very low temperatures 
(T/J ~ 10~ 2 ) enabled us to be as close as possible to 
the asymptotic regime where universal scaling laws are 
expected to hold. 

Figure n shows the specific heat of CS2CUCI4 in zero 
magnetic field. The magnetic contribution, C mag to the 
total specific heat, Ctot, was obtained by subtracting the 
phonon contribution C ph = 13599 (T/@ D ) 3 Jmol"^- 1 
(using a Debye temperature 0£> = 126 K). Furthermore, 
the contribution of the nuclear specific heat to the to- 
tal specific heat at very low temperatures has been ac- 
counted for in the a naly sis of all the raw data shown 
in the following 0, lll| . The two prominent features 
present in C(T) are the broad maximum related to the 
cross-over from the paramagnetic to a short-range spin 
correlated state and the A-like peak. The latter is the 
signature of the entrance into the 3D magnetically or- 
dered state (Tjv = 0.595 K), where the magnetic struc- 
ture is a spiral in the (6, c) plane |3j . The overall shape 
of C(T) above Tjv is already captured quantitatively by 
including only the strongest term in the Hamiltonian, 
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FIG. 1: Specific heat of CS2CUCI4 in zero magnetic field. 
From the total specific heat, Ctot (open symbols) the phonon 
contribution C v h has been subtracted to obtain the magnetic 
part Cmag (filled symbols). The solid and dashed lines rep- 
resent the calculated C(T) (see text). Inset: Temperature 
dependence of the magnetic entropy, S m a g (T). 
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FIG. 2: Magnetic specific heat, C mag vs T of CS2CUCI4 close 
to the critical field. Inset: Cmag/T vs T in a semi-logarithmic 
plot. The data recorded at B = 8.44 T reveal a small jump 
at T c — 76 mK, indicating the vicinity of the critical field. 



the dominant coupling J (solid line) Exact theo- 

retical predictions are not available for the full Hamil- 
tonian in CS2CUCI4, but including the next order term, 
the frustrated 2D zig-zag coupling J' = J/3, and using 
a high-temperature series expansion technique |l3j. give 
similar behavior (dashed line) although with a slightly 
lower specific heat at the position of the broad maximum. 
It may be possible that including the other small terms 
like the DM interaction (D a /J ss 5%) could improve the 
agreement. The magnetic entropy S mag (T) shows that 
S = i?.ln2 for spin- 1/2 is almost recovered at 10 K (in- 
set to Fig. ^ which corresponds roughly to the observed 
band width of the spin excitations at base temperature 
0. 

The ordering temperature and the position of the 
broad maximum hardly change for small fields. However, 
the transition temperature to the spiral ordered state, 
which can be regarded as a cone-like structure £| , varies 
very strongly above 8 T (Fig. Instead of T/v, we la- 
bel the transition temperature in this field range T c , as 
in eq. ^ m order to follow the nomenclature used in the 
theoretical description. The A-like anomaly in C rraQg (T) 
is gradually suppressed in its height and its position is 
pushed to lower temperatures with increasing field. An 
extraordinary change occurs as the field is increased from 
8.4 T to 8.44 T (inset to Fig. HJ. Upon this tiny field 
change (AB/B < 0.5%), T c is reduced by almost a factor 
of two (T c = 76 mK at B = 8.44 T), and T c has shifted 
downwards by almost one order of magnitude compared 
to the zero-field value. No further signatures of the tran- 



sition can be resolved in our data upon approaching the 
critical field B c ~ 8.5 T. For B > B c the ordering of 
the transverse component of the magnetic moment com- 
pletely disappears since the spin system enters a field- 
induced ferromagnetic (FM) state 

A field-induced gap in the magnon excitation spec- 
trum was observed in the FM state by neutron scattering 
measurements Q. Its field dependence was given to bc 
A = gfi B (B - B c ), with B c ~ 8.44 T, g = 2.18, and fi B 
the Bohr magneton. For the interpretation of the phase 
transition below B c as a BEC of magnons it is crucial 
that the gap closes at B c . To provide a compelling evi- 
dence for this fact we re-examined the phase diagram £| 
above B c with our thermodynamic measurements. The 
results are presented in Fig. |3| 

The magnon dispersion along the a-direction is small 
due to a weak interlayer spin coupling. Thus, for tem- 
peratures well above a characteristic energy scale E* « 
50 mK, the actual magnon dispersion is of 2D charac- 
ter. However, for T < E* a smooth cross-over to a 3D 
character is expected. Assuming first a 2D quadratic 
magnon dispersion, the leading contribution to the tem- 
perature dependence of the specific heat is given by 
C m ag — exp(— A/T)/T, provided that T < A. As shown 
in Fig. this behavior fits well the experimental data 
above 0.3 K. The obtained field dependence of A(B) is 
discussed below. The deviation from a straight line of 
the 9 T data below 0.3 K in Fig. [3] might indicate the 
cross-over from 2D to 3D magnons. This notion is sup- 
ported by the low-temperature data plotted as C mag VT 
vs 1/T in the inset to Fig- EI This presentation is used 
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FIG. 3: Semi-logarithmic plot of CmagT vs 1/T of Cs 2 CuCl 4 
for fields above B c . In this representation the slope of the 
data (solid lines) yields the value of the gap A present in the 
magnon excitation spectrum. The data shown in the inset 
were obtained at B — 8.6 T and are plotted as C ma gVT vs 
1/T. 



since a 3D dispersion relation yields as leading term in 
the specific heat C mag ~ exp(— A/T)/v^T. The straight 
line below w 0.15 K indicates that this model describes 
the data appropriately. However, in this cross-over re- 
gion, the 2D model describes the data equally well and 
the determined value of the gap is (within the error bar) 
the same as the one deduced from the C mag T vs 1/T 
plot. 

The (T,B) phase diagram of CS2CUCI4 obtained from 
our specific heat experiments is presented in Fig. The 
field dependence of the magnetic transition temperature 
is in very good agreement with T C {B) obtained from neu- 
tron data up to 8 T ^4|. The specific heat data, how- 
ever, revealed that T c starts to decrease strongly above 
8 T and T c — ► for B — ► B c . Fitting the power law de- 
pendence T C (B) oc (B c - B) 1 ^ to the data for B > 8 T 
with the assumption of B c — 8.50 T yields an exponent 
<p = 1.52(10). We want to stress that the value of <\> is 
very sensitive to the chosen value of B c . An exponent 
4> = 1.44(10) is obtained if B c = 8.51 is used. The solid 
line plotted in the inset to Fig. ^represents the result of 
the theoretical analysis described below. Above B c the 
fully spin-polarized FM state is created and the gap A 
opens in the spin excitation spectrum. The dashed line 
represents a linear fit to the data for B < 10 T (main part 
of Fig. 0. This yields B c = 8.3(10) T and g = 2.31(15). 
The relatively large errors are due to the uncertainties in 
the fit. 

To treat the observed phase transition slightly below 
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FIG. 4: (T, B) phase diagram of Cs 2 CuCl 4 for B \\ a. The 
ordering temperature T c decreases as the magnetic field ap- 
proaches the critical field B c — 8.51 T. Above B c a field- 
polarized ferromagnetic (FM) state is entered and the gap 
A in the spin excitation spectrum opens. Inset: The experi- 
mental T C (B) data points are well described by the calculated 
phase boundary of the BEC of magnons (solid line) . 



the saturation field B c as a BEC of magnons 0, Q, @ 
we used the hard-core boson representation for spin- 1/2 
operators S^, Sf in the original Hamiltonian H. Because 
the DM interaction (D = 0.053(5)J 4]) changes sign 
between even and odd magnetic layers, which are stacked 
along the a-direction, two types of bosons, a,, and bj are 
introduced for the two types of layers [jja, [l6| ■ The hard- 
core boson constraint was satisfied by adding to H an 
infinite on-site repulsion, U — > 00, between bosons given 
by 



H. 



(a) 



H, 



(b) 



U H h PPi h i- (2) 



The interlayer coupling J" — 0.045(5) J mixes a and 
b boson modes and results in two bare magnon excitation 
branches A and B. Their dispersion relations are given 
by i 



E^ B — J q T BignD q ^Dl + (j;') 2 - E , 

J q = J cos q x +2J' cos (fe/2) cos(q y /2) 
D q = 2Dsm(q x /2)cos(qy/2) , 
j" = j"cos(q z /2) . 



with 



(3) 

(4) 
(5) 
(6) 
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Here J' = 0.34(3) J £| and the (/-values are restricted to 
< q x < 2tt, ^ q y < 4tt, and ^ q z < 2tt. 

The degenerate minima E A - = E B - = are at Qi = 

Qi Qi 

(n + 6i, 0, 0) for branch A and at Q2 = (tt — 82, 2tt, 0) for 
branch B. Without loosing precision we can use <5i ~ 
82 — 8 — 2 arcsin ( J'/2 J). The bilinear part of H now 
reads 

Hut = y, [K - w>) A t A i + ( E * - wj) B i B i] > ( ? ) 

with A q = a q a q + (i q b q , B q = a q b q - (3 q a q , a 2 + f3 2 = 
1, and the bare chemical potential fiQ — g[iB (B c — B). 
The saturation field B c — W/(gfiB), with W being the 
magnon bandwidth, was calculated to be B c = 8.51 T 
assuming g = 2.20 fl6lll7|. 

The interaction given by eq. |21 describes the scattering 
of A and B magnons. Near the quantum critical point, 
(B c — B) <C B c and at low temperature, the average den- 
sity of magnons n A = n B = n is low, n ~ (1 — B/B c ). 
The magnon scattering can be treated in the ladder ap- 
proximation 01 1 neglecting interference between a and 
b channels. In this approximation, the problem reduces 
to solving the Bethe-Salpeter equation in each channel. 

This results in the renormalized scattering amplitudes 
rW(gi,^;g 3 ,g 4 ) for i = a,b. Here cf 3 ,g 4 and qi,q 2 are 
magnon momenta before and after scattering, respec- 
tively, and qi + q\ = <f 3 + q^. The total energy of scat- 
tered magnons was set to zero. This limit is compati- 
ble with our main goal to describe the phase transition 
near B c when approaching the phase boundary T c (B) 
from higher temperatures. At T — > T c , only the magnon 
states at q ~ Qi,2 are occupied and the magnon spectrum 
renormalization near the minima is important. 

With given r( a ) and T^ b \ the complete set of two- 
particle scattering amplitudes was then obtained by mul- 
tiplying r'°)'' b J by products of four a q and [3 q coeffi- 
cients. For instance, a scattering process (A^ 3 ,B^ 4 ) — > 
(Ag\ , Bq~ 2 ) in the channel a is described by the amplitude 
a qi (3 q2 a q3 f3 q j: {a) {qi, q 2 ; q 3 ,qi)- 

The renormalization of low-energy magnons was found 
by treating the magnon scattering effects in the Hartree- 
Fock approximation: 

H™ F = 2TnY{A + q A q + B+B q ) 
1 

+2T'a(3nJ2 {^+B q + B+ A q ) , (8) 
where Oq i = ag 2 = a 2 and 0q 1 — Pq 2 = P 2 ■ Taking 



into account that a 2 ^> 1 , we keep here only the leading 
contributions to energy parameters T and V: 

= a 4 r^ (g 2 , g 2; q 2 , q 2 ) , r' ~ 2r/a 2 (9) 

and we obtained the estimate T ~ 0.85J. According to 
eq. |S1 the chemical potential of magnons is renormalized 
/Xo — * Me// = Mo — 2rrt, and the low-energy magnons 
are hybridized due to the second term in eq. |H1 This 
term shifts the bottom of the magnon band slightly down 
and leads to a weak mass enhancement of low-energy 
magnons. Both effects are proportional to n 2 and we 
omit them since n <C 1 near T c for (B c — B) <C B c . 

For a given B < B c and with decreasing temperature 
the magnon BEC occurs when the effective chemical po- 
tential fJ-eff vanishes 6]. Then T C (B) is determined by 

gUB (B c -B) = 2Tn (T c ) . (10) 

Here n(T) = with f B (E q ) being the Bose 

distribution function taken at fJ. e ff — and E q = E q 
or E q — E q . This means that for T < T c the magnon 

condensate develops simultaneously at q = Qi,2- It is 
worth emphasizing that at fJ-eff ~ > the distribution 
function fsiE) diverges as T/E for E — > 0. Therefore, 
the low energy 3D-magnon spectrum, E < E* , mainly 
contributes and drives the BEC transition. The phase 
boundary can be calculated using eq. 1101 It gives a very 
good description of the experimental data near B c (see 
inset to Fig. QJ, but deviates strongly at lower fields, 
i.e., for B c — B > 0.5 T. This indicates that the mean- 
field description of the magnon BEC is only applicable 
in the close vicinity of B c . The calculated boundary is 
well described by eq-Qwith a critical exponent 4>th — 1-5 
close to the predicted value 0bec = 3/2 characteristic for 
3D quadratic dispersion of low-energy magnons 0, |(| . 

We have presented experimental evidence that in 
CS2CUCI4 the field dependence of the critical temper- 
ature T C (B) oc (B c - B) 1 ^ close to the critical field 
B c = 8.51 T is well described with cf> ~ 1.5. This is in 
very good agreement with the exponent expected in the 
mean-field approximation. Together with the observed 
opening of a spin gap above B c these findings support 
the notion of a Bose-Einstein condensation of magnons 
in Cs 2 CuCl 4 . 
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